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dxy- like pairings in organic superconductors k-(BEDT-TTF)2X 
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Pairing symmetry in «:-(BEDT-TTF)2X is re-examined by applying fluctuation exchange method 
to the Hubbard model on the original four band lattice for the first time. Our study is motivated 
by the discrepancy between a recent thermal conductivity /tunneling spectroscopy experiments and 
the previous theories based on the effective single or two band model. In contrast to the previous 
theories, we find that d^2_y2-like pairing dominates over d^y in agreement with recent experiments. 
Our study provides a general lesson that even if the Fermi surface seems to be reproduced by an 
effective model, the pairing symmetry can be different from that of the original system. 

PACS numbers; 74.20.Rp, 74.70.Kn, 74.20.Mn 

Unconventional mechanisms of superconductivity is 
one of the most challenging issues in solid state physics. 
Not only that there is a possibility of accomplishing high 
Tc as in the cuprates, but also the variety of pairing sym- 
metries provide rich physics. 

Recently, it has become increasingly clear that organic 
materials can provide various stages for unconventional 
pairings. In particular, k-(BEDT-TTF)2X is of special 
interest in that superconductivity lies near the antiferro- 
magnetic insulating phase like in the cuprates, ^ sug- 
gesting that spin fluctuations are playing important role 
in the occurrence of superconductivity. 

In order to clarify the mechanism of superconductivity, 
it is crucial to identify the pairing symmetry. Theoreti- 
cally, previous studies on the Hubbard model with fluc- 
tuation exchange (FLEX) approximation (^-^ or quan- 
tum Monte Carlo calculation |^ have suggested that su- 
perconductivity in K-(BEDT-TTF)2X is similar to the 
d-wave superconductivity in the cuprates. However, a 
recent thermal conductivity as well as tunneling spec- 
troscopy measurements casts doubt on such a view, 
which is exactly the motivation of the present study. 

Let us first briefly summarize the background of the 
issue. As shown in Fig.|l|(a), k-(BEDT-TTF)2X is essen- 
tially a four band system in the sense that there are four 
BEDT-TTF molecules per unit cell. The lattice structure 
is dimerized in that the hopping integral in the 5-direction 
alternates as tbi, ib2, i&i, • • where > \tb2\- The band 
filling is n = 1.5, where n is defined as n =(the number 
of electrons/ the number of molecules). The upper two 
bands, which we will refer to as a and (3 bands hereafter, 
cross the Fermi level, so the Fermi surface consists of two piQ i ^j^^ j^^^j^^ structure of k-(BEDT-TTF)2X. (a) the 

portions (see Fig|2|). In the limit of large hi, namely, original four band lattice, and (b) the effective two band lat- 
when the dimerization is strong, each dimer can be con- tice. bl, fo2, ■ ■ ■ stand for the hopping integrals t^i, t62, ■ ■ ■. A 
sidered as a single lattice site, so the system reduces to unit cell consists of the hatched molecules (or sites). 



a two band model shown in Fig.0(b) with n = 1, where 
the band filling is now defined as n=(the number of elec- 
trons/the number of sites). In other words, the en- 
ergy gap between the upper two and the lower two bands 
becomes large when the dimerization is strong, so that 
the lower two bands, which do not cross the Fermi level, 
can be neglected. In this limit, the effective hopping inte- 
grals tb and tc are given as tb — — tb2/2, = (— tp + iq)/2, 
and tc' — {—tp' + tq')/2 in terms of the original hopping 
integrals, ||] which gives |tb/tc| 0.8. The system fur- 
ther reduces to a single band model when tc — tc' . 

Assuming that the dimerization is sufficiently strong in 
K-(BEDT-TTF)2X, the two band |] or the single band 
Hubbard model [^[js) has been adopted in the previous 
theoretical studies. There, it has turned out that an- 
tiferromagnetic spin fluctuations lead to an anisotropic 
superconductivity with a gap function of the form shown 
in Fig.^(a), which we will refer to as da;y-like pairing here- 
after from the direction of the nodes. This gap function 
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FIG. 2. (a) di:y-like and (b) d^'2_y2-\Tike gap functions on 
the Fermi surfaces. Solid (dashed) lines represent positive 
(negative) values. The inset in (a) shows the Fermi surface 
and the daj^-like gap function in the unfolded Brillouin zone, 
which can be obtained by moving the right and the left por- 
tions of the P band Fermi surface in the folded Brillouin zone 
by (— 27r,0) and (-|-27r,0), respectively, and rotating it by 45 
degrees. Arrows show the position of the nodes of the gap 
function. Note that the points where the a and the /3 band 
Fermi surfaces (almost) touch in (a) are not the nodes. There, 
the gap function abruptly jumps from a negative to a positive 
value (or vise versa). This can be understood by looking at 
the corresponding points in the unfolded Brillouin zone. 
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is similar to that of the d-wave superconductivity in the 
high Tc cuprates, as can be seen more clearly by unfolding 
the Brillouin zone (inset of Fig||(a)). 

Experimentally, NMR experiments |^^|-|l2| have indeed 
supported the spin-fluctuation scenario, but those exper- 
iments do not give direct information on the position of 
the nodes in the gap function. A millimeter-wave trans- 
mission experiment has suggested a da;j,-like gap function, 
|p3| but different interpretations on this experiment have 
recently been proposed, 15| so the situation has not 



been settled. 

Quite recently, tunneling spectroscopy [Q and thermal 
conductivity measurements |^ have been performed for 
(BEDT-TTF)2Cu(NCS)2 in an attempt to determine the 
anisotropy in the gap function. Surprisingly, both of the 
experiments have suggested that the gap function is more 
like Fig||(b) (referred to as d2;2_j,2-like hereafter) than 
(a) , casting doubt on the spin-fluctuation scenario based 
on the single band or the two band approximation. 

In the present study, we apply FLEX to the Hubbard 
model, H = Y.^,^ E^=T,i ^yklcj^, + H_c) + [/ ^ . n,|n,X, 
on the original Jour hand lattice (Fig.n^(a)) for the first 
time. Surprisingly enough, our results show that d^2_j,2- 
like pairing dominates dxy when dimerization is not so 
strong, as in k-(BEDT-TTF)2Cu(NCS)2. The present 
study provides a general lesson that even if the Fermi sur- 
face seems to be well reproduced by an effective model, 
its pairing symmetry can be quite different from that of 
the original system. 
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FIG. 3. Contour plots (for the lowest Matsubara fre- 
quency) of (a)|G''p, (b) <^^._^., (c) <f>ly (v = a, /3) and (d)x 
with tbi = -1.0, = -0.5, tf, = -0.45, tq = 0.15, [7 = 6, 
and r = 0.01. The thick lines in (b) and (c) represent the 
nodes of the gap function. 



FLEX is a kind of self-consistent random phase approx- 
imation, which is known to be appropriate for treating 
strong spin fluctuations. ||l^-[l^ In the four band version 
of FLEX, [|9| the Green's function G, the susceptibility 
X, the self-energy S, and the superconducting gap func- 
tion (/) all become 4x4 matrices, e.g., G/m(k, where 
Z, m specify the four sites in a unit cell. The site- indexed 
matrices for Green's function and the gap functions can 
be converted into band(denoted as a, /3, • • •)-indexed ones 
with a unitary transformation. As for the spin suscepti- 
bility, we diagonalize the 4x4 matrix x&v ^'^'^ concentrate 
on the largest eigenvalue, denoted as x- We take 64 x 64 
/c-point meshes and up to 4096 Matsubara frequencies. 
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FIG. 4. A^2_j^2 and Xxy plotted as functions of temperature 
for the same parameter values as in Fig.^ 



The superconducting gap function is obtained as the 
eigenfunction of the linearized Eliashberg equation. In 
order to solve the Eliashberg equation, we use the power 
method, where da;2_y2-like or da;y-like gap functions has 
been assumed as a trial function. Consequently, we have 
obtained d^^-y^ and da;y-like gap functions 0^2_j,2 and 
= Of; ■ ■ respectively, together with the corre- 
sponding eigenvalues A2;2_j,2 and Xxy 

We have modeled k-(BEDT-TTF)2Cu(NCS)2 by tak- 
ing till — —1, tb2 — —0.5, tp — —0.45, and tq = 0.15 in 
accordance with a band calculation result given in ref. 
pO| . (We have neglected the small difference between tq 
and tqi or between tp and tpt.) In Fig.^ (a), we show 
contour plots for |G"^p {v = a,/3), where the ridges rep- 
resent the Fermi surfaces. In (b) and (c), we show con- 
tour plots for (f>'^2_y2 and (p^y, respectively. Superposing 
the gap functions on the Fermi surfaces, it can be seen 



that I 



and ' 



bxy essentially have the functional forms 
given in Fig.^ (b) and (a), respectively. 

The corresponding eigenvalues Aj.2_j,2 and X^y are 
plotted against temperature in Fig.^. The tendency that 
Xx2-y2 — Xxy increases upon lowering the temperature 
clearly demonstrates that dx^-y^-^^^e pairing dominates 
over dxy 

The reason why dx2-y2-\ike superconductivity dom- 
inates over dxy can be found in the spin structure. 
Namely, the spin susceptibility x, shown in Fig.^(d), 
peaks around Q ~ (±0.57r, ±0.67r), which is a conse- 
quence of partial nesting between the a and the /3 band 
portions of the Fermi surface (FigJ^(a)). Thus, the inter- 
band pair scattering processes dominate over intraband 
ones, resulting in a gap function that has different signs 
between the two portions of the Fermi surfaces, but does 
not change sign within each portion. 
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FIG. 5. (a) The nesting vector Q bridges the a and the P 
band portions of the Fermi surface when the dimerization is 
not so strong. (b)The nesting vector lies within the a band 
portion of the Fermi surface for the single band and the two 
band models. The same applies for the four band model when 
the dimerization is strong. d^2_y2- and d^^y-like gap functions 
are shown in (a) and (b), respectively. 




FIG. 6. Contour plots of jCp and x for the single band 
model with it = 0.83, ^ 1, U ^ 6 and T = 0.03. Folded 
Brillouin zone is adopted for clear comparison with FigJ3 



The present FLEX results for the four band model is 
indeed different from those for the corresponding single 
band model. In Fig.^, we show |G|^ and x for the single 
band Hubbard model with n = 1, tb = 0.83, tc = tc' = 1, 
and U = 6, where the ratio tb/tc is taken to be the same 
as [-tb2/2]/[{-tp + tq)/2] with 42 = -0.5, tp = -0.45 
and tq — 0.15. (We have adopted the folded Brillouin 
zone in order to make clear comparison with the results 
for the four band model.) Although the Fermi surfaces 
look similar to that of the four band model (Fig.||(a)), the 
spin susceptibility peaks around (0, ±0.257r) in contrast 
to that of the four band model. Thus the pair scattering 
mainly occurs within the a band portion of the Fermi 
surface (Fig.^(b)), and consequently dxy-^i^e supercon- 
ductivity dominates in the single band model. A general 
lesson here is that even if the Fermi surfaces of an effec- 
tive model 'look similar' to that of the original system. 
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FIG. 7. (a)A contour plot of x for the same parameter val- 
ues as in Fig.^ except that tti — —1.6. (b)A^2_j^2 and \xy 
plotted against \tbi\ at T — 0.01. 

the spin structure, and consequently the superconducting 
gap function, can be quite different. 

When the dimerization is strong, the four band model 
should tend to the single band model. In fact, if we take 
for example thi = — 1.6 and leave the other parameters 
unchanged, x, shown in Figj7|(a), now has a structure 
similar to the single band result Fig.^(c). Consequently, 
(ij,2^2-like superconductivity gives way to d^y as seen in 
Fig.[7|(b), where we have plotted Aj.2_j,2 and Xxy against 
\tbi\ at T = 0.01. 



As we have seen, the competition between d^^ 



and 



(ia;y-like pairings is sensitive to the parameter values. 
This means that the two types of superconductivity can 
possibly be nearly degenerate for certain combinations of 
the anion X, temperature, pressure, and/or magnetic 
field. In such a case, dj,2_y2 -\-idxy pairing, a pairing with 
broken time reversal symmetry, is highly probable since 
the gap is fully open in this channel. Investigations on 
such a possibihty would be interesting both theoretically 
and experimentally. 

To summarize, we have examined the Hubbard model 
on the original four band lattice of k-(BEDT-TTF)2X. 
We have found that dx^-y2-\\ke superconductivity dom- 
inates over dxy when the dimerization is not so strong, 
as in K-(BEDT-TTF)2Cu(NCS)2. Our finding is con- 
sistent with recent thermal conductivity and tunneling 
spectroscopy measurements. 
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The spin structure of k-(BEDT-TTF)2X is known to be 
commensurate with the lattice in and near the insulat- 
ing phase This is because the localized spin picture 
is valid in the insulating phase. However, when the sys- 
tem is metallic and the Fermi surface is well defined (as 
is expected near the superconducting phase), the spin 
structure should be affected by the nesting properties of 
the Fermi surface and in general can be incommensurate. 
For the parameter values adopted in the present study, 
d^2_y2 clearly dominates over d^y. However, since there 
are some ambiguities in the estimation of the param- 
eter values, we believe that there is still a possibility 
that d^2_y2- and d^y-like pairings closely compete with 
each other (under certain conditions) even for the case of 
X=Cu(NCS)2. It is interesting to look, from this view- 
point, at a recent experiment suggesting the existence 
of a phase transition within the superconducting phase 
under magnetic field [J. Singleton, J. A. Symington, M.- 
S. Nam, A. Ardavan, M. Kurmoo, and P. Day, J. Phys. 
Cond. Matt. 12, L641 (2000).] 



